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TIME-AVERAGES OF FAST OSCILLATORY SYSTEMS IN 
THREE-DIMENSIONAL GEOPHYSICAL FLUID DYNAMICS AND 
ELECTROMAGNETIC EFFECTS 

BIN CHENG AND ALEX MAHALOV 


Abstract. Time-averages are common observables in analysis of experimental data and nu¬ 
merical simulations of physical systems. We will investigate, from the angle of partial dif¬ 
ferential equation analysis, some oscillatory geophysical fluid dynamics in three dimensions: 
Navier-Stokes equations in a fast rotating, spherical shell, and Magnetohydrodynamics subject 
to strong Coriolis and Lorentz forces. Upon averaging their oscillatory solutions in time, in¬ 
teresting patterns such as zonal flows can emerge. More rigorously, we will prove that, when 
the restoring forces are strong enough, time-averaged solutions stay close to the null spaces of 
the wave operators, whereas the solutions themselves can be arbitrarily far away from these 
subspaces. 
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1. Introduction 

In many geophysical fluid dynamical systems, solutions exhibit fast oscillatory behaviors due 
to strong energy-preserving, restoring mechanisms — a typical example being the Coriolis force 
in fast rotating planets and stars. Time integration averages out the oscillatory part of the 
solution, which leads to emerging of interesting patterns that are relevant in a longer time scale. 

A straightforward framework is introduced in [7] for proving that the time-average of the 
solution stays close to the null space of the large, skew-self-adjoint operator in the partial 
differential equation (PDE) system. A particular application of this framework can be found in 
[6] for two-dimensional (2D) Euler equations on a fast rotating sphere. 

In this article, we study two PDE systems in three-dimensional (3D) spatial domains that 
are important models of fast oscillatory, geophysical fluid dynamics. 

The first one, (2.2) - (2.3), governs viscous, barotropic fluids confined within a fast rotating, 
spherical shell that models the global atmospheric circulation on Earth and other planets. 
On the boundary, the velocity held either satisfies conditions in terms of shear stress or is 
simply fixed. We prove in Theorem 2.1 that, with additional spatial-averaging in the radial 
direction, time-averages of the solution are 0(e) close to zonal flows (i.e. motions in the east-west 
direction). Here, e denotes the Rossby number, a dimensionless parameter measures the ratio 
between typical magnitudes of inertia and Coriolis force. This theoretical result is consistent 
with many numerical studies and observations. For a partial list of computational results, we 
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mention [13, 23] for 3D models, [33, 8, 25, 15, 30, 12] for 2D models, and references therein. 
Note that many of these computations attempt to simulate turbulent flows with sufficiently 
high resolutions. Zonal structures in these numerical results are either directly noticeable by 
naked eyes or after some time-averaging procedures. 

On the other hand, we have observed zonal flow patterns (e.g. bands and jets) on giant 
planets for hundreds of years, which has attracted considerable interests recently thanks to 
spacecraft missions and the launch of the Hubble Space Telescope (e.g. [14], [27]). In [24], the 
banded structure is directly observable in a composite view of the Jovian atmosphere captured 
by the Cassini spacecraft. There are also observational data in the oceans on Earth showing 
persistent zonal flow patterns (e.g. [28, 29, 22]). 

It is worth mentioning that time-averaging and the more general time-filtering are commonly 
used for denoising of observational and computational data. Indeed, such post-processing is 
necessary for the emergence of zonal flow patterns in some of the above literature, e.g. [23, Fig 
4,6 and 9]. We also point out that the zonal flow pattern would not arise in a model without the 
meridional variation of the Coriolis parameter. Such variation is due to the non-flat geometry of 
the spatial domain, which is why we use the entire spherical shell as the domain in this article. 
Most analytical work in literature adopts the /3-plane approximation, focusing on a narrow strip 
near a fixed latitude, which essentially is a linear approximation of the spherical case. Fourier 
series then become applicable in the /3-plane approximation but is not so in the whole spherical 
shell that is studied in this article. 

Mathematical studies of deterministic and stochastic 3D rotating Navier-Stokes equations 
including resonances were done in [2, 3, 4, 11] with uniformly large rotation. Recently, [17] 
proves some interesting /3-plane effects using the randomly forced quasi-geostrophic equation. 

The second PDF system (3.1) governs rotating Magnetohydrodynamics (MHD) in the whole 
three-dimensional Euclidean space subject to two strong restoring forces: Coriolis force and 
Lorentz force. They induce the magneticstrophic waves, also known as rotating Alfven waves 
([10]). We refer to [21] for ionospheric applications. We prove in Theorem 3.1 that time-averages 
of the solution vanish at order of fractional powers of e when measured in L* norms (s > 6). 
This result suggests that there is dispersion in the time-averages, although we do not impose 
any spatial decay on the initial data like in the classical dispersive wave theory. 

The rest of this article is organized as following. The formulations and main results are 
introduced in Sections 2, 3. Then, in Section 4, we apply the barotropic averaging (2.5) on 
the 3D Navier-Stokes equations (2.2) and reveal the close connection to the 2D Navier-Stokes 
equations on a sphere. In Section 5, we prove the main Theorem 2.1 using the time-averaging 
tools devloped in [6, 7]. In Section 6, we study the MHD system (3.1) and prove Theorem 3.1 
by using Sobolev-type inequalities. Finally, in Section 7 the Appendices, we give a geometric 
proof of Proposition 2.3 regarding the Navier boundary conditions and also prove energy and 
enstrophy estimates for the Navier-Stokes equations in a shell subject to Navier boundary 
conditions with A > 0 and g = 0. 
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2. Navier-Stokes equations in a rotating shell: formulation and main resnlts 

Let (x, y, z) denote the usual Cartesian coordinates and let spherical coordinates (r, 9, 4>) 
denote the radius, colatitude (i.e. inclination from the positive half of the 2 ; axis) and longitude 
respectively. The spatial domain is a thin, spherical shell 

H := {{x,y,z) \ \/x'^ + y"^ + z'^ e (1 - 5,1 J)} 

= {(r,6',(/>) I r G (1 - (5,1 -h (5),6» € [O,7r],0 G [0,27r]}. 
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( 2 . 1 ) 



In short, we can write 0 = (1 — (5,1 + 5) x with denoting the unit sphere. 

Let er,eg,e^ denote the locally orthogonal unit vectors along the increasing directions of 
r,9,(j) respectively — and they are orientated according to the right hand rule, i.e. (e^ x e^) -e^ = 
1. Similarly define ex,ey,ez in terms of the Cartesian coordinate system. The unknown is 
velocity field u. The Coriolis force is given by 


1 ^ 

Fcoriolis = “U X 6^ COS 9 = -U X 6^ 
e s 

where e, called the Rossby number, equals the ratio of the spatial domain’s rotating period over 
the inertial time scale (usually 0.01 ~ 0.1 for the Earth). Note we have adopted such geophysical 
version of the Coriolis force that differs from the laboratory version, x e^. In other words, 
we neglect the radial component of the velocity and also neglect the radial component of the 
Coriolis force. See [35] for detailed justification. 

Let q denote the pressure and constant the viscosity. The sum of other external forces 
is denoted by Fext- Then, the incompressible Navier-Stokes equations under the Coriolis force 
reads ([1, 9, 26]) 

+ VuU + Vg = -u X er cos 9 + /xAu + Fext, 


div u = 0, 

subject to the Navier boundary conditions which consist two parts, 

u • = 0, (zero-flux) 

Sn -I- Au = g, (given shear stress) 


- tan 


( 2 . 2 ) 

(2.3a) 

(2.3b) 


an 


with scalar A = X{t, x) > 0 and vector g = g(t, x) given. (The physical significance of having 
positive A is shown in Proposition 7.1 and its proof.) Here, n denotes the outward normal at 
clH and subscript “tan” indicates the tangential component, e.g. vtan := v — (v-n)n. The stress 
tensor S is defined as 


where 


S := Vu + (Vu)T 


/ dxUi 

dyUi 


Vu := dxU2 

dyU2 

dzU2 

\dxU3 

dyU3 

dzU3/ 


Throughout this article, vectors are treated as 3 x 1 matrices so that for vector fields u, u', u^'. 


(Vu)u' = u'Vu, 


(2.4a) 


((Vu)Tu') • u" = (u"-Vu) • u'. (2.4b) 

One can also impose the Dirichlet boundary condition, = g' with g' given. This appa- 
rantly includes the non-slip boundary condition = 0. 

Before stating the main result, some definitions are in order. First, the Sobolev norm for 
a scalar or vector function / is defined as 

ll/llL2(n) := ^ I/p• 

Second, we use the so-called “baratropic averaging” to reduce the 3D velocity field u to a 2D 
field u that is tangent to It turns out that certain weight in the integral is convenient. From 
a physical perspective, the flux of u going through a side of the area element sin 9 d9d(f> equals 
the (unweighted) radial average of the momentum flux through the corresponding vertical cross- 
section of the volume element r'^ sin 9 dddtpdr. To this end, for velocity field u S L^(D), define 
its horizontal component 

U/j := u — (u • 6^)6^ 
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and define its baratropic averaging 


_ 1 

ruhdr. (2.5) 

Ji-s 

Next, define U : ^ L^{S'^) as the zonal-mean projector that projects horizontal 

zonal 

velocity fields onto the subspace of zonal flows, 

Uu{e, [ u{9,(l)) ■ e^d(i] e^. 

zonal yZTT Jq J 

Lastly, let C denote some universal constant and we add subscript(s) to it, e.g. Cfc, to 
emphasize its dependence on another parameter. 

We now state the main result subject to homogenous boundary condition g = 0. We will 
skip the case of nonhomogenous boundary condition g 7 ^ 0 because of the intimate connection 
between g and the external forcing Fext that is discussed in Subsection 2.2 where such connection 
is explained also in a more physical context. 


Theorem 2.1 (Homogeneous boundary conditions). Consider 3D Navier-Stokes equations (2.2) 
in a spherieal shell 12 defined in (2.1), subject to the Navier boundary conditions (2.3) with 
g = 0. Let pL < lj2 and 5 < 1/2. Define 

Mo ■= ;^^I|uo||l2(q), 
that indicates the averaged size of uq . 

Then, for any weak solution n G L°°([0, 00); L^(n))nL^([0,00); Ll^(n)), its barotropic average 
u as defined in (2.5) satisfies, for a < —4, 


(1 


n 


zonal 



i7“(S2) 


< e [(1 + C/rT)Mo + CaTM^ + M^^t ], 


( 2 . 6 ) 


where Mgxt 



n )/<, r 


ext 


/i'“+2(S2) 


and the norm can be defined using spherical 


harmonics (c.f. Definition 5.1 and relation (5.5). Note for negative a, the H°‘ norm dampens 
high wave number modes). The constant Cn depends solely on a and is otherwise independent 
ofe, 6, fi, Mo, T. 


Since operator (1 — H ) effectively extracts the non-zonal component of a velocity field, 

zonal 

T 

estimate ( 2 . 6 ) confirms that Jq u is 0 (e) close to zonal flows. 

Combining this theorem with the energy estimate (7.11), we can obtain via interpolation 
that, for a G [—4,1), 




ndt 


H“(S2) 


< [(1 + CfiT)Mo + CaTM^ + Mext] eV 


a .—b 


where a, b are positive numbers depending on a. 

We remark that the possible negative a values used in the estimates above suggest that 
zonal flow patterns are associated with smaller wave numbers i.e. larger spatial scales, since the 
high wave number modes are damped in the definition of norms for negative a. 

The above results in 3D are nontrivial extension from the 2D case studied in [ 6 ] which is 
centered around the Euler equations on a fast rotating unit sphere 

clju-I-VuU-|-V(? =-u-*-, divu = 0, 
e 

where denotes the 7 r /2 counterclockwise rotation of the associated vector on §^. 
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(2.7) 
















For comparison, the main theorem for the 2D system (2.7) is stated as following with some 
minor notational changes. 


Theorem 2.2 ([6]). Consider the incompressible Euler equation (2.7) on with initial data 
uq G 77^(S^) for k > 3. Let Mq := ||uo||/^fc. Then, there exists a function /(•) : [—1,1] M 
depending on u, so that 


(u - V-^/(z)) dt 






( 2 . 8 ) 


for any given T G [0, r*/Mo] where constant T* depends on k but is independent of e and uq. 
In spherical coordinates, 

V-^f{z) = -/'(sin 6 ')sin 6 »e 0 , 
which represents longitude-independent zonal flows. 


We finally remark that analysis of 3D Navier-Stokes equations and its variations in the 
geophysical context, including the existence of solutions and the low Rossby number limit, has 
seen substantial progress in recent years, e.g. [19, 18, 20, 2, 3, 4, 5, 11], just to name a few. There 
are also results regarding Navier-Stokes equations on thin 3D domains, e.g. [32, 16], without 
the Coriolis effect. The boundary conditions in the literature are either periodic, whole space, 
non-slip = 0 or some variations of the Navier type (2.3). In the next two subsections, we 
further discuss (2.3) and its variations as seen in literature. 


2.1. Geometry of the Navier boundary conditions. The following proposition is regarding 
a general domain D. 


Proposition 2.3. For general smooth domain Q C M^, let n be the outward normal at a point 
of dQ. Suppose 

u-^lao = 0- 

Then, the Sn term in the Navier boundary conditions (2.3) satisfies 


[S'n]tan 


dn 


[n 'Vu]^^^ — n -Vn 
(curlu) X n — 2 u -Vn, 


(2.9) 


where the u-Vn is well defined (intrinsically) on dQ due to u-n|g^ = 0. We can further rewrite 
it using 

2 

- u-Vfi = '^{u-ei)Kiei, (2-10) 

i=l 

with 6 * 1,62 being a pair of orthonormal, principal directions of the surface dLl and k* = — (e* • 
Vn) • e) being the corresponding principal curvature. 

As a consequence, for the case of spherical shell domain defined in ( 2 . 1 ), 


[5*11] tan 


= zb 

r=ldi(5 


duh 

dr 



= (curl u) X n =p 


2 ufe 

r 


( 2 . 11 a) 

( 2 . 11 b) 


The proof is postponed to the Appendices. Also, consult [34] for more details. 

In literature, the Navier boundary conditions are also referred to as “stree-free” or “slip” 
boundary condition. It should be however distinguished from the so-called “free” boundary 
condition (which is confusingly referred to as “slip” boundary condition in some cases). 


u • n|gQ = 0, (curlu) X ^1^^^ = 0. 
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The Navier boundary conditions should also be distinguished from the Neumann type boundary 
condition, 

duhi 


u 


^\dn 


= 0 


Qj. Ian 

as used in e.g. [18]. By (2.11) of Proposition 2.3, they only correspond to special cases of 
the Navier boundary conditions (2.3) when one lets g = 0 and makes specific choices for the 
parameter A. We are not aware of any physical explanation for these choices of A. A more 
serious issue is that such specific choices always involve A < 0 on the inner boundary r = 1 — 6. 
However, by Proposition 7.1 and its proof, the constraint A > 0 on the boundary everywhere is 
necessary for the dissipation of energy. Also it is physically invalid to argue that the above two 
boundary conditions are the small-curvature approximations of the Navier boundary conditions 
(2.3), because the principal curvatures of Sfl are of 0(1) in global circulation models for which 
the radius of Earth is rescaled to near the unit. 

Interested reader can further consult [16], in particular the top part of page 1085, and refer¬ 
ences therein. 

2.2. Physical considerations of external forcing and non-homogeneous boundary 
conditions. In the main Theorem 2.1, the external force Fgxt affects the estimate only via 
(1 — n ) Fext which is its non-zonal component averaged in time and r. This external 

force is intimately connected to non-homogeneous boundary conditions which are studied in 
e.g. the context of planetary boundary layer (PBL). Mathematically speaking, if u satishes the 
Navier-Stokes equations (2.2) with nonhomogeneous boundary conditions 


u • n 

S^n + An 

J tan 


dfi 

= g / 0, 

an 


and if one can find some velocity field v, regardless of the dynamics, that is only subject to the 
boundary conditions 

/ -*1 
V • n 

on 

( 2 . 12 ) 


Svfi + Av 


J tan 


an 

= g, 


an 


Then, the new unknown u := u — v will satisfy (2.2) with homogeneous boundary conditions 




S^n + An 


J tan 


= 0 , 


an 


which is then covered by the main Theorem 2.1. The new external force term in the u system 
apparently contains information of the original boundary data g. 

There are indeed infinitely many ways to construct v satisfying (2.12). For example, it suffices 
to hnd vector fields a(0, 4>), b(0, cj)) that are both tangent to so that 

V = r^a -I- b 

satisfies (2.12). The v • = 0 part is apparently valid. For the second condition of (2.12), 

we rewrite it using (2.11a), 




r=ldi5 


= g 


r=l±(5 


Substitute v = r^a -|- b and rearrange 


-1 


(l -I- (5 -|- (1 -|- 5)^A)a -|- (^ ^ ^ -|- A)b — g 
( — l + (5 + (l — 5)^A)a-|- (^^—- -|- A)b = g 


r=lH-5 


r=l—S 
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With A > 0 and <5 G (0,1/2), we always have 

(1 + 5 + (1 + < 5 ) 2 a ) > I - 1 + 5 + (1 - < 5 ) 2 a | 


and 


■1-5 


+ A) > 


-1 

TT5 


+ A|. 


Therefore, the coefficient matrix of the above linear system is diagonally dominant. Thus, we 
can perform Gaussian elimination (while pretending a, b to be scalar unknowns) and express 


a, b as linear combinations of g 


r=l+(5 


g 


which are both tangent to 

r=l—5 


3. Magnetohydrodynamics in formulation and main results 

Consider the domain to be in which a uniform, imposed magnetic field e^ resides and 
a fast rotating (about e^), conducting fluid moves subject to the predominantly large Coriolis 
force and Lorentz force. The fluid is homogeneous, incompressible and un-magnetizable. Then, 
upon some scaling arguments, one can reduce the full Navier-Stokes and Maxwell’s equations to 
the following MHD system [10, §3.8] for the unknowns: velocity field u and induced magnetic 
field b (so that the total magnetic field is given by + eb), 


dtu + u -Vu + Vq 
u X e. 


(curl b) X (e^ + eb) 


divu = 0; 


dth = 


curl [u X (e^ + eb)] 


divb = 0. 


(3.1a) 


(3.1b) 


Here, e denotes the MHD Rossby number as well as the ratio of the induced magnetic field 
over imposed magnetic field; q denotes the pressure. Note in (3.1a) the Coriolis force and Lorentz 
force are of the same scale which is 0(l/e) times the inertia. This is a reasonable scaling since 
the ratio of these two forces is often close to 1 in many geophysical and astrophysical applications 
([10]). For simplicity, we have set both the kinetic viscosity and magnetic viscosity to be zero. 

Let Mo := ||(uo, bo)||/^fe(R3) for k > 5/2. By the standard energy method, we know 


ll(u, b)||j^fe(]g3) < CkMo 

for positive times t < C^/Mq. 


(3.2) 


Theorem 3.1. 

T < Ck/Mo, 


and 


Consider any classical solution to (3.1) satisfying (3.2). Then, for any positive 

<CkMo[{T + MoT^)e]K (3.3) 


/ 


udt 


3,00 


f 


hdt 


Wfc-4.s(M3) 


<Ck,sMo[{T + MoT^)e 


(3.4) 


with 6 < s < oo. 


Therefore, time-averages of the solution vanish at order of fractional powers of e when mea¬ 
sured in norms (s > 6). This result suggests there is dispersion in the time-averages, although 
we do not impose any spatial decay on the initial data like in the classical dispersive wave theory. 
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4. Barotropic averaging of the Navier-Stokes equations 

Recall the definition of barotropic averaging (2.5), 

. 1+5 


- W w 

^ TTT / dr, 

2o Ji-s 


(4.1) 


where 

u/i := u — (u • e^) e,.. 

Also define the barotropic average for a scalar /, 

1 


(4.2) 


We hrst remove the pressure term Vg in (2.2) using the Helmholtz-Leray decomposition. 
Dehne X to be the space of incompressible velocity helds subject to zero-flux boundary condition, 

X := closure of 


|u”" G I divu'“ = 0 , = o} 


By using testing functions, we see that 
X = closure of 


u G ‘^^(11) f u -V/ = 0 for any / G 

Jq 

= |u G L^{Vt) J u -V/ = 0 for any / G ff^(fl)| 


(4.3) 


Dehne V as the L^-orthogonal projection onto X so that, for any u, u' G L^(D), 

V^u = pu G X, 


/ (u - Vu) ■ {Vu') = 0. 

Jq 


(4.4a) 

(4.4b) 


In fact, V is the classical Leray projection subject to zero-hux boundary condition. Then, dehne 

Q:=I-V. 

Now pick any scalar / G By orthogonality of V, Q in (4.4b), we have 


/ V{Vf )-V{Vf )= [ lP(V/)-V/ 

Jq Jq 

which is zero due to P(V/) G X satisfying (4.3). In other words, 

V{Vf) = 0. 

By this property, we apply V on the hrst equation of (2.2), cancel the Vg term and reformulate 
it into, 

dtu + P(Vuu) = -P(u X cos 6) + ^VAu. (4.5) 

Note that, for generic div-free velocity held u satisfying the Navier boundary conditions (2.3), 
the term An is no longer subject to the zero-hux boundary condition u • = 0 whereas the 

image of V always satishes the zero-hux boundary condition. Thus, VAw and An differ by a 
div-free, potential how — the gradient of the so-called Stokes pressure. 





4.1. Barotropic averaging of Helmholtz-Leray projection. Similar to (4.3), we define 
X/i = |u/j G J Ufe-Vc/= 0 for any 5 G 

and then define projections Vh and Qh := I — Vh for “horizontal” velocity field G L^(§^) so 
that, analogous to (4.4) 

'Ph^h = Vh^h G Xh, (4.6a) 

[ Qh^h ■ i'PhW) = 0. (4.6b) 

J §2 

Here and below, subscript h following an operator indicates the operator acts on scalar or 
vector fields defined on S^. For example, A/j denotes the LaplaceBeltrami operator on S^. Their 
properties are discussed in the appendices of [7] . 

Now, we give the relation between V and Vh- 

Lemma 4.1. For any vector field u = Vn + Qu G L^(n), 

Vvi = VhVi, Qu = Qhu. 

In the proof, we will repeatedly use the following basic facts that relate the differential 
operators in Q to those in 
For vector u = wBr + Uh, 


d 

divu = r~‘^ — {r‘^w) + r” 
dr 

curlu = r“^(curl/jU/i)ej. 

d 


div/jU/i, 


(4.7a) 


+ r ^{VhW - —{ruh)) X 


for scalar /, 


Vf = ^fer + r-^Vhf, 
dr 


(4.7b) 


Note that the relation of curl and curl/i in polar coordinates (with 6 being the colatitude) is 
due to the following formulations that roughly resemble the Cartesian-coordinate form, 

for u = wer + U 0 eg + u^e^, 

1 


curl u = 


r sin0 
1 

+ - 




+ 


J_A 

r \^sin0 d(p 
1 


l(rne) - ^(u.) 


curl,,u/j = ( ^(u^sin^) - -^ug] 

smtl \dd o<P J 


Proof of Lemma f.l. Apply barotropic averaging (4.1) to u = Vvl + Qu and get u = Vvl + Qu. 
Since by definition we also have u = VhVi + QhU, h suffices to prove 

Qu = Q/iU. (4.8) 

Also, since is dense in Lp‘{Q) and barotropic averaging is apparently bounded from LfiQ) 

to T^(§^), we only consider u G 77^ (Q) so that u • n is defined on dPt. 

By elliptic PDF theory, we have 

A/ = divu in Q 

V f ■ n = u ■ fi in dLl 
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Qu = V/ where 


(4.9) 






Similar equations hold for (c.f. [7, (2.1), (2.2)]), 


QhVih = V/jA^ ^ dwh^h (4.10) 

Here, is defined using spherical harmonics, and maps between scalar functions of zero mean 
— note /g 2 div/jUft = 0 by Stokes’ lemma and = 0. 

Let u = wSr + u/j. Use (4.7) to reformulate (4.9) as, 

-2 ^ / 2 ^ p\ , -2, 


or or 


= r 


-2 


dr 


(r^w) + r div/iU/i, 


<9 f 
or 

Multiply the first equation with and integrate it in r, 
d i+<5 ^ i+<5 


in H, 
on do,. 


(4.11) 


r—f 

dr-’ 


+ / Ahf dr = r w 
l-<5 Ji_5 


rl+S 

+ / rdivhUhdr. 
i-<5 Ji_5 


Then, apply the second equation of (4.11) to cancel out the boundary terms, 

^l+<5 rl+S 

/ Ahf dr = r divhUh dr. 

Ji-5 Ji-5 

Since u € 77^ (H), we can exchange integrals and derivatives above, and invoke definitions of 
barotropic averaging in (4.1), (4.2) to obtain 

A/i7=divftU, i.e., V hj = Qh^, (4.12) 

where Qh follows (4.10). 

On the other hand, apply barotropic averaging (4.1) on the first equation of (4.7b) with the 
same / as in (4.9) to obtain 

U+5 


Qu = V/ = 


^Jl-S 


Vhfdr = Vhf. 


Combine it with (4.12) on S^, we prove Lemma 4.1. 


□ 


4.2. Dynamics of barotropic averages on EP. We now apply barotropic averaging (2.5) on 
the 3D Navier-Stokes equations (4.5) with the help of Lemma 4.1 and identities (4.7). 

Lemma 4.2. The solution to (4.5) subject to the Navier boundary conditions (2.3) with g = 0 
satisfies 

1 


dtu + P/iVuU = -'Ph(u X e^ cos^) + nVhAu + Fext 

subject to div/jU = 0. Here, u = we^ + Uh so that u = vifi. 

Furthermore, the viscosity term from above equals 


VhAu = VhAhV + 2Vhr 

Proof. First, integrate divu = 0 in r and invoke the first identity of (4.7a) 

1+5 ri+<5 

1-5 


(4.13) 


(4.14) 


r'^w 


rl+S 

+ / rdiYh{uh)=d. 

Ji-5 


By the zero-flux boundary condition u • = 0, the first term vanishes, and therefore we 

prove the incompressibility condition 

div/jU = 0. (4-15) 

For the Coriolis term, V{u x e^cos^), Lemma 4.1 implies 


V{u X e^ cos 6) = P/i(u x e,. cos 6) = P/i(u x e^ cos 6). 
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Then, upon barotropic averaging and invoking Lemma 4.1, the 3D Navier-Stokes (4.5) is 
transformed into (4.13) subject to div/^u = 0. 

Now we show (4.14). By divu = 0, we have, 

Au = —curl curl u. (4-16) 

For the RHS, first apply the second identity of (4.7a) to get 

U+5 


_ 1 D+'’ Q 

curlu = / {S/hW - —{ruh)) x dr 

=4(i 


VhW dr — ru/i 


Then, substitute u by curlu and correspondingly substitute w = (curlu) ■ = r ^curl/jU/^, 


— curl curl u 
-1+5 


1 ( f 

'25 VA-5 


Vhr ^curl/iU/i dr — rProj/j(curlu) 


1+5 


1-5 


X e,, 


= X Vhcuilhr-'^uh + — Proj;j(curl u) 

= :/ + //. 

Thus, we transform the viscous term in (4.13) into 


1+5 

1-5 


X Sr 


P/iAu = —Rftcurl curlu = Vhl + VhH- 
For the Vhl term, apply Vh on the identity 

AfeU/j = Or X V/iCurl/iU/j + divh^h 
and then use Vh^h = 0 to rewrite 


(4.17) 


Vhl = Vh^hr ^u/i. 


(4.18) 


For the VhH term, invoke the second identity of (4.7a) and the zero-flux boundary condition 


u • fi|g^ = 0 to obtain 


so that. 


1+5 _ Q 1+5 

Proj/j(curlu) =-r — (ru/^) x 
1-5 or 


1-5 




1+5 

1 - 5 ' 


By (2.11), the Navier boundary conditions (2.3) with g = 0 imply = r ^u/j. Therefore, 

1+5 


II = luj 


1 d 

= 2r ^^u. 
1-5 or 


Combine this with (4.17), (4.18) to prove (4.14). 


□ 


5. Proof of main theorem for Navier-Stokes equations 

In this section, we following the framework in [7] to prove Theorem 2.1. 
First, define 

■='Ph{u X Or COS 6) 

and rewrite (4.13) as 


(5.1) 


ChU = e 


dtu + VhV^u+ 


-2^Vhr ^drVi- liVh^hr - Fext 
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Then, take the time-averages of each term and exchange time integration and 

Ch 


j u — e [u(T) — Uq -|- -|- ^2 -|- ^3 ~ Fext 


(5.2) 


where 


: = 


"P/iVaUdf 


A 2 — 2/i f 

Jo 


Vhr~^drudt 


"L 


^3 := -Ai / VhJ^hr "^Vihdt 


We will then estimate every term in the RHS of (5.2) in terms of norms. Note that 

there are many equivalent definitions of Sobolev norms on a manifold through the literature 
(e.g. [31]), all of which are independent of coordinate systems. One definiation of ||/||//fc(s 2 ) is 

f/ Ei=o(-^)^'/l for ^ > 0- Then, by the Poincare’s inequality, this definition is equivalent 


adopt the following definition 


. In this article, all relavent scalar fields are of zero-mean, so that we 


for scalar / with / / = 0 , 
J §2 


dehne ||/||hRS 2 ) := 




(5.3) 


with integer /c > 0 . 

Consequently, for a vector field u on with Hodge Decomposition 


u = V/i<h-|-with / = / 

J §2 


= / T = 0, 


we define its norm, among other equivalent versions, as 


|n|lH''(S 2 ) • Y^II‘^’IEfc+l(§2) + l|T|Ei,+i^g2y 


(5.4) 


Note that, here and below, we always impose zero-mean on and T. 

With the help of spherical harmonics, we extend (5.3) and (5.4) to for any a G M. 


Definition 5.1. Let {L)™} for ^ = 0,1,.. . and m = —I, ..., —1,0,1,..., / be the set of spherical 
harmonics forming an orthonormal basis of such that 

Ahvr = -i{i + i)Yr. 

Let u = be any div-free velocity field in L^(§^) with 

00 I 

l=l rn=—l 

Then, for any real number a, 


|u||//a(g 2 ) 


= 114' 


00 I 


Ih“+i(s2) := a E E 

\ i=i 


l=l rn=—l 


This definition allows us to easily adapt the proof of Theorem 4.1 of [ 6 ] and reach the next 
lemma (whose proof is skipped). 
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Lemma 5.2. Let a G M. For any horizontal vector field on subject to div/jU/j = 0, 

11(1- n )u,,||j^a(§2) < ||£/jU,j||j^a + 2(S2). 

zonal ^ ^ ' 

Note that (1 — 11 ) effectively extracts the non-zonal component of a velocity field. 

zonal 

Proof of Theorem 2.1. Let Ll stand for the three-dimensional shell domain defined in (2.1) for 
the rest of the proof. We will also use without references the integrating-by-parts formulas on 
which can be found in e.g. [7, (A.19)-(A.22)]. 

Under Definition 5.1, it is straightforward to verify that if“(§^) and 77““(§^) are dual spaces 
with respect to the L^(S^) inner product, namely. 


|u||i/-a(s 2 ) — rn 


max 


(U,u')i 2 (§ 2 ) 


uVO ||u||jya(g 2 ) 

Then, by Lemma 5.2 and (4.15), it suffices to estimate 


(5.5) 


/ u,u,' 
Jo 


L 2 (§ 2 ) 


for smooth, testing vector field that is tangent to S^. Since the definition (5.1) implies 
div/jT/j = 0 , we can further impose div^Uh = 0 so that for any u/^" tangent to S^, 


{'Ph'^h ) W )l 2('§2) — (W ,U/i )i^2(^§2y 


(5.6) 


By (5.2), it suffices to make the following estimates. (Recall definition Mq := ||uo||j;^ 2 (q)/\/^.) 
• Estimate of u(T) and ug. 


2S(u{T),Uh' 

\ /L 2 (S 2 ) 

rl+S n 

= / ru(r) • Uh 

Jl-S 

< l|u(T)||L2(f2)||r"^u,,'||i2(t2) 

= l|u(T)||L2(Q)||u/j'||i2(§2)A/^ 

Since ||u(r)||j;^ 2 (f 2 ) < ||uo||j;^ 2 (q) = v^Mg by Proposition 7.1, we obtain 

I|u(^)IIl 2 (§ 2 ) < ^ 0 , and similarly ||ug||i 2 (s 2 ) < Mg 
• Estimate of Ai. By (5.6) 


2 ^V 

rl+S r- 

- / r [V • (u (g) u)] • u,,' 

Jl-5 


= (V • (u(g) u),r ^u/i' 
= (u ^ u,V{r~^Uh)^ 


L 2 (Q) 


1,2 (n) 


by = 0 


< l|u|lL2(f^)|V(r u/,')Il-(q) 

<C'l|u||i2(j^)|Uft'|wl,°c(§2) 

< C'/3||u||i2(Q)||Uft'||jy2+/3('§2) 


for /3 > 0 


. by Sobolev imbedding. 
13 




Since ||u(r)|| 2 , 2 (f 7 ) < I|uo||l 2 (q) = V^Mq by Proposition 7.1, we use the duality relation (5.5) 
to obtain 

||'PhVuU||j;^-2-/3('g2) < 

and integrating in time gives 
for /3 > 0 . 

• Estimate of ^ 2 - By (5.6) 


ll^l|liT-2-/5(S2) < CpTMo^ 

rl+S /■ 

Jl-5 7s2 


(9rU • Uh' 


= (drU,r ‘^Uh 


Then, integrate in time to get, 


L 2 (n) 

< l|Vu||i2(f2)||r“2u;,'||^2(t2) 

< C'v^||Vu||j;^2(q)||u/i'||^2(§2). 


^ 2 , U/i' 


L 2 (§ 2 ) 


^ , f 

< I|U/i'||l 2 (§ 2 ) ||Vu||i 2 (f 2 ) 


< 




Since by the enstrophy estimate (7.11) 


l|Uft'||L2(S2) 


llVul 


Il 2 (o)- 


iff 


l|Vu||2,^^^ < ^WnoWHl+ CfiT) /fi 


we obain 




II^2||l2(§2) < C -\/+ CfiT) Mq 

Estimate of A^. By (5.6) and Green’s identity on (c.f. [7, (A.22)]) 

= 26 (r-‘^Uh, AhUh 


r-l+<5 


'L 2 (S 2 ) 

rl+S r 

■■ / / r"^u/j • AhUh 

Jl-S J§2 


,-2, 


= {r ^Ufe,r ^AhUh' 

< C'I|u||l2(o)||A/iU/i'||j;^2(§2)V^. 

Since ||u(T)||j;^ 2 (o) < ||uo||l 2 (o) = V^Mq by Proposition 7.1, we obtain 

||'ff’/iA/ir“ 2 u||j^- 2 (§ 2 ) < CMo, i.e. ||A 3 ||j^- 2 (§ 2 ) < C^TMq 


□ 
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6. Fast Rotating MHD model 

Let V once again denote the Leray-Helmholtz projection. In other words, for any vector field 

U€ L2(m3)^ 

u = Vvl + VQu 

where 

Pu = —curl A”^curlu, Qu = A“^divu. 

By Fourier transform, we also have 


= - 


X X u) 


Qu(0 = 


■ u 


|iCP' 

Dehne a skew-self-adjoint operator acting on the velocity field u and magnetic field b. 


( 6 . 1 ) 

( 6 . 2 ) 


Af 


P(u X 02 -|- (curlb) X e^) 
curl (u X e^) 


Then, the system (3.1) can be reformulated as 


df 


u 




T’(uVu-bVb) 
u Vb — b Vu 


= -Af 
e 


(6.3) 


(6.4) 


where we used identity 2(curl b) x b — 2b -Vb = V|bp to transform (3.1a) and identity curl (u x 
b) = u( divb) — b( divu) -I- b - Vu — u - Vb to transform (3.1b). 

6.1. Kernel of the large operator It is an elementary calculation to verify that 

divu = 0 curl (u x e^) = SzU. 

Combine it with (6.1) and the fact divu = divcurlb = 0 to transform (6.3) into 


Af 


—curl A ^^^(u-|-curlb) 


d, 


u 


Thus, by b = Vh = —curl A ^curlb, this implies 




—curl A ^dzVi + d^b 
dzVL 


(6.5) 


Therefore, 


(u,b) £ ker..Sf u = u(x,y) and b = h{x,y) 


When restricted to space, kerAf = {(0,0)}. 

6.2. Control L°° norm using dz derivatives. Given initial data {k > 5/2), we can apply 
the standard energy method to obtain 0(1) estimates for the U’ norms {2 < p < oo) of the 
solution and its first spatial derivatives in a finite time interval. Upon time-averaging, such 

0(1) estimates give rise to 0(e) estimates for the norms {2 < p < oo) of Jq Af ( ^ j. By 

(6.5), the 0(e) estimates also work for (S^u, d^b) in terms of norms^. 

What estimates can be obtained for (u, b), the time-average of the solution itself? Because 
of the special role of dz derivatives, we state and prove the following inequality regarding function 
/ defined in M^. 

ll/lli°°(Ri) ^ C\\f\\L^(R^)\\f'\\L^{R^) (6-6) 

Note that once the L°° norm is estimated, standard interpolation techniques can help control 
the rest of the norms {2 < p < oo). 


^For dzh, the range of p is reduced to 6 < p < oo due to the negative derivative in the first line of (6.5) 
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Proof of (6.6). For any positive p, we estimate the norms of /(^) over frequencies lower and 
higher than p respectively. The Holder’s inequality is applied in both cases. 


I'P ^ r rP ii:rrP 

/ |/(0|de< / l/(OP / 1 

'—p ^J—p ^J—p 


< C’II/IIl2(ki)P2 . 


[ \m\df< 


1 

2 


h^\>p 

^Ai\>p 



< C’ / L2(r1)P 2 




Therefore, 

ii/iiloo(mi)< [ \m\df 

^ C'||/||i2(Kl)p2 + C'||/'||i2(Rl)P 2. 

Optimizing the RHS over p G (0,oo), we prove (6.6). □ 


I 2 
Iff” 


We are ready to state and prove the following lemma 
Lemma 6.1. Given function g(x, y, z), 

1 

II5||lj°(h™(k2)) < C\\g\\f^^^^3^\\dzg\ 

Consequently, for and m' > 0 and k > 1, 

1 1 

llff II — ^lll^ll H'^'^fl'll ■ 

Proof. Consider d^dy g{x,y, z) with a + ft < m. For any numbers a < b, we estimate 

/ \\dxd^9\\l2/z 

J CL 


b — a 


b — a 


[d^d^g{x, y, z)] ‘^dxdydz 


< 


< 


\\dxd^9{x,y,z)\\l^^.^i^dxdy 
^ J ^2 ^ {d^dfd^gf dz ^ dxdy 


... by (6.6) 

<C ( {pfd^gYdzdxdy 

- J R3 - - J Jlf^ 

... by Holder’s inequality in 

< C'llfl'll//"*(R3)l|52fl'llfl'"*(K3) 


[dzdfd^gf'dzdxdy 


Because a, b are arbitary, this implies 

ll'9“^y5'llLoo(i2^) < C'||p||j^m(R3)||0^fi(||jym(K3). 

Summing up over all derivatives with 0 < a + /3 < m, we complete the proof of the hrst 
inequality and the second one follows from the Sobolev inequalities. □ 
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6.3. Estimates on time-averages of (u, b). 


Proof of Theorem 3.1. By taking on (6.4) and then taking the H 
the estimate 

< CMo ( + Mo 1 e. 


fc-i 


!>(;) 


jjk-i 


norms, we have 
(6.7) 


f 


Using the second component of (6.5), we obtain from above that 

< CMo -|- Mq^ e, 

and together with Lemma 6.1, we prove (3.3) 

Similarly, substracting the second component of (6.5) from the curl of the Hrst component 
of (6.5), we obtain from (6.7) that 

rT 


f 


52Curlb||j:^;=-2(K3) < CMq — + Mo e 


and consequently, by Lemma 6.1, 


11 curl 


b||^fc-4,ooC]g3') < C*Mo ( — + Mo 62. 


( 6 . 8 ) 


To “remove” the curl operator from (6.8), we use b = Vh = —curl A ^curlb (—A) zb = 
curl (—A)“ 2 curlb and the fact that curl (—A )“2 is a bounded mapping on any 1T™'’^(M^) space 
with 1 < p < oo (by properties of Fourier multipliers) to obtain 


-A)l f 
Jo 


T r-T 

b||l4/m,p(]]j3) < C*||curl / b|| 54 /m,p 

0 Jo 


Then, apply the Hardy-Littlewood-Sobolev fractional integration theorem || 5 r|| 3p < C\\ (—A) a^H^^p 

2 l_l 

to the LHS and the interpolative Holder’s inequality \\g\\Lp < Ils'lli 2 lbllL°°^ RHS to arrive 

at, with 2 < p < 3, 

rT 

b|| 3p 

/o 


< 


(711 curl f 
Jo 


bll ^ 

^\\\Ym,2 


11 curl 


b ^ 


Finally, plug in (3.2) and (6.8), we complete the proof of (3.4) by setting m = k — 4, 1/p = 
1/3 + 1/s. □ 


7. Appendices 

7.1. Geometric proof of Proposition 2.3. 


Proof. Throughout this proof, let n denote the outward normal at dQ and let r denote a typical 
tangent vector at dH. 

By identities (2.4), we have. 


(7.1) 


J ((Vu)n) • T = (n-Vu) • T, 

I ((Vu)'''n) • r = (r-Vu) • ft. 

Next, the assumption u • =0 implies r • V(u • = 0 so that by treating r • V as a 

directional derivative and using the product rule, we have at dH, 

(r • Vu) ■ n = -{f-Vfi) 


u. 
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Combine it with (7.1) to obtain, at a general smooth boundary dQ with u • = 0 


((Vu)n)-r= n-Vu • r, 

((Vu)'''n) • r = —{r-'Vn) ■ u. 

Recall the definition S = Vu + (Vu)t and apply the above identities to obtain, 


{Sn) ■ T = 


n Vu 


T — {t -Vn) ■ u. 


(7.2) 


(7.3) 


One can also write the above identity using the vorticity. In fact, combine identities (7.2) with 
(curlu) X n = (Vu — (Vu)''')n to have 


n • Vu 


r + (r • Vn) • u. 


(7.4) 


((curlu) X n) ■ T = 

Subtract it from (7.3) to arrive at 

{Sn) ■ T = ((curl u) X n) • r — 2(r • Vn) • u 
Now, regarding the -f-'S/n term, it is associated with the shape operator^, 

S^{f) := -T-Vn 

which is a linear mapping in any given tangent plane of dVl. Then, the symmetric bilinear form 

^{ri,T 2 ) := y{Ti) • T 2 = -(ti - Vn) • T 2 , (7.5) 

defined for any two tangent vectors Ti,r 2 is the the second fundamental form^. Its symmetry 
can be shown straightforward e.g. by choosing the surface as a level set of scalar function g, 

so that n = which we will skip. Combining such symmetry with (7.3), (7.4) and the 

assumption u • = 0, we prove (2.9). 

To prove (2.10), we recall that the (orthonormal) principal directions^ are the two (orthonor¬ 
mal) eigenvectors of the shape operator 


y{ei) =-Si-Vn = Kiii, i = l,2. 


(7.6) 


where eigenvalue Hi denotes the principal curvature associated with Cj. They effectively diago¬ 
nalize i.e. by definition (7.5), 


{^i i^j) — 


0, i + j- 


Combine (7.6), (2.9) to obtain (2.10). 

Finally, for the special case of II being a spherical shell, any pair of orthonormal tangent 
vectors can serve as the principal directions, and therefore (2.11) easily follows from the fact 
that n = iOj. at r = 1 ± (5 and from (7.6) so that 

Ki = —{Si -Vn) • Si = =F- at r = 1 ± d. 

r 


□ 


^http;//mathworld. wolfram.com/ShapeOperator.html 

^http: / / mathworld.wolfram.com/SecondFundamentalForm.html 

'^https://en. wikipedia.org/wiki/PrincipaLcurvature 
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7.2. Energy and enstrophy estimates in a thin shell with Navier boundary condi¬ 
tions. Let norm || • || stand for the L‘^{Q) norm and (•, •) for the L‘^{Q) inner product. 

Proposition 7.1. Consider (2.2) subject to the Navier boundary conditions (2.3) with A > 0 
and g = 0. Then, the energy ||u||(t) is decreasing with time. 

Proof. Take the L^(fl) inner product of u and the first equation of (2.2), noting the Coriolis 
term is perpendicular to u, 

/i(Au, u) = (dfU, u) + (u • Vu, u) + (Vg, u) 

= (^u -Vlul^ + u -Vq) 

2 " " 

where the last step is due to the Divergence Theorem, zero-flux boundary condition u-n|g^ = 0 
and div u = 0. 

Now, it is useful to derive a version of the Green’s formula adapted to the Navier boundary 
conditions. First, use divu = 0 to rewrite An = div S' so that 


(An, u) = (div S, u) = —(S, Vu)-|- f (Sn) ■ u. 

Jdn 

The second term is non-positive. In fact, at dD, the Navier boundary conditions (2.3) imply 

that u is perpendicular to n while (Sn -|- An) is parellel to n for g = 0. Therefore, with A > 0, 

/ (Sn-|-An) • u = 0 so / (Sn) • u < 0. (7.8) 

Jdn Jdn 

For the (S, Vu) term, use the definition of inner-products between matrices and the fact S^ = S 
to obtain, (S, Vu) = (S''', (Vu)''') = (S, (Vu)'''). Therefore, 

(S,Vu) = i(S,Vu + (Vu)T) = i||Sf. 

So, combine the above 3 equations to arrive at 

(Au,u) < -^||Sf. 

Together with (7.7), it implies 

dt||u||^ < -/r||S|p. (7.9) 

The proof is complete. □ 


To obtain some estimates on the total enstrophy ||curlu||, it suffices to estimate ||Vu|p. 
Simply taking the time integral of the above inequality (7.9) will however not yield estimate on 
||Vu|pdt because S lacks the information on the anti-symmetric part of Vu (which actually 
conincides with curlu). The remedy is to employ another version of Green’s formula 


(An, u) = — ||Vu|p-|- / (n-Vu)-u 

Jdn 

= -\\Vuf+ {n-Vuh)-Uh, 

Jan 


(7.10) 


where the zero-flux boundary condition u • = 0 was also used. For the boundary term 

above, apply the first equation of (2.11) and the fact that n = Te^ at the boundaries r = 1 ± A 
to obtain, 


n-Vu/i 


= [-S'njtan + {n-er) — . 
an r 
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Then 


/ (n-Vu/i) • u/j = / 

JdQ. Jan 


< 


tan T ■ 6 r) ) ' U/j 

r J 

J ^[5n] + (n • • u 

[ (n • e^)- • u • • • by (7.8). 

Jan f 


Effectively, there are no more derivatives of u in the boundary integral (indeed this formula 
works for general domain). Thus, by applying the divergence theorem to the right side above 
to obtain, 

/ (n-Vu,,) • u,, < / div(er|upr"^) 

Jan Jn 

< C (||Vu||||u|| + ||u|p) 

where we also used the Holder’s inequality. Substitute it into (7.7), (7.10) to obtain 


-dt\\uf < ||Vu||2 

1 , 


C'(||Vu||||u|| + llul 


</r(^-^||Vuf+ C||uf 

with a different constant C. Together with the decrease of energy ||u IP < l|uo|P due to (7.9), 
it implies 


/ 


IVufdt < iiuolpf- + cr). 
V/r 


(7.11) 
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